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Coherent states of the two dimensional harmonic oscillator are constructed 
as superpositions of energy and angular momentum eigenstates. It is shown 
that these states are Gaussian wave-packets moving along a classical trajec- 
tory, with a well defined elliptical polarization. They are coherent correlated 
states with respect to the usual cartesian position and momentum operators. 
A set of creation and annihilation operators is defined in polar coordinates, 
and it is shown that these same states are precisely coherent states with 
respect to such operators. 
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1 Introduction 



Coherent states were introduced by Glauber jT] in quantum optics and have 
a very wide range of application in quantum physics. They can be defined in 
several equivalent ways: as Gaussian wave packets evolving without spread- 
ing, and also as close analogues of classical states satisfying the minimum 
dispersion relation allowed by the Heisenberg uncertainty principle. One 
of the various generalizations of coherent states are the coherent correlated 
states introduced by Dodonov et a/.|2], that satisfy the minimum Robertson- 
Schrodinger uncertainty relation. 

A coherent state can be expressed as an infinite superposition of the 
number states of a harmonic oscillator. The generalization to two or more 
modes can be easily achieved by coupling several linear harmonic oscillator. 
This is the usual approach used in this kind of problems; for instance, Arickx 
et al. jHj have studied the evolution of Gaussian wave packets in two and 
three dimensions in cartesian coordinates. 

However, polar coordinates are the natural coordinates for a two-dimensional 
oscillator with a single frequency, permitting to separate directly the angu- 
lar momentum. Accordingly, it would be advantageous to have a definition 
of coherent states that takes into account the rotational symmetry of the 
problem. This is particularly relevant due the recent interest in the angular 
momentum of light, associated to electromagnetic fields with axial symmetry 
such as Bessel and Laguerre beams (see e. g. Padgett et al. [4J). 

In this article, we propose a new kind of coherent states that are superpo- 
sitions of states with well defined energy and angular momentum, and that 
can be directly related to axially symmetric fields. The analysis of these 
states is presented in Section 2, where it is also shown that they are coher- 
ent correlated states with respect to cartesian position and momentum, and 
correspond to two independent elliptical polarizations, rotating clockwise or 
counterclockwise. In Section 3, a set of creation and annihilation operators in 
polar coordinates are defined, and it is shown that the states are truly coher- 
ent with respect to these operators because they have minimum dispersion. 
Some numerical examples are presented. 
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2 Coherent correlated states 

Consider a two dimensional harmonic oscillator with mass M and frequency 
uj. Choosing M, to -1 , and {h/muj) 1 / 2 as units of mass, time and length 
respectively, the time dependent Schrodinger equation in cylindrical coordi- 
nates takes the form: 

<i»M,») = \{~ ;|r| - ~ + r*)nr,M- (2.1) 

The eigenstates of energy and angular momentum are \n, ±Z) (we take I > 0). 
In coordinates representation, these states are given by the standard text- 
book solutions (see, e. g., [5J): 

(rA\n,±l) = (-j^) 1/2 e ±u *e-r 2 /\%(r 2 ), (2.2) 

where L l n {r 2 ) are Laguerre polynomials. The energy of these states is E n j = 
2n + l + l. 

In analogy with the usual definition of coherent states, let us now define 
the following states: 



\a,(3) ± = N(a,[3) £ {^±Jl) 1/2 ^ a n p l \n,±l), (2.3) 

n,l=0 ^" 



where N{a 1 (3) is a normalization constant. Clearly, the evolution of these 
states is given by 

e~ iH \a, P)± = \ae~ 2it , pe~ u ) ± . (2.4) 
The orthogonality condition determines N(a,/3), that is: 

00 (n + IV 

±(a,P\a,fl± = \N(a } f3)\ 2 £ ^J-±\ a n^ 1 = 1. (2.5) 

n,l=0 n \ L ) 

Using the formula 

y ^±H x n - n (26) 



it follows that 



with the condition |a| < 1. 
It can also be seen that 



±(a',{3'\a,{3) ± 



a 



a: 



x 



exp{ 



1 — aa*' 

,|2 pp, 



2(1 -\a\ 2 ) 2(1-|«'| 2 ) l-aa* 



}, (2-8) 



and 



: (af,/3'\a,l3) ± = 



a 



a' 



x exp | — 



1 — aa*' 

\P\ 2 \PT 



2(1 - |o;| 2 ) 2(1 - |«'|2) J • (2 ' 9) 

Let us now look for the explicit form of the states \a,(3)± in coordinates 
representation. Since 



1 00 1 

n,i=0 



7T 



(2.10) 



and using the generating function formula: 



^ a n L^(x) = (1 - a)- , - 1 e x °/( a - 1 ), 

n=0 



(2.11) 



it follows that 

(r,0|a,/3)± 

Accordingly: 



-(l+o)r 2 /2(l- 



iV ± (a,/?) 
7rV2(i-a; 

^(r,0) = (r,0|a,/3) ± 



(2.12) 



1 /l + a r 1 + a r 2 fie** 

-^VT^ exp i i^ T + T^ 



}. (2.13) 



a 2(1 — | a 

Now, writing the exponent in this last equation in the form: 

- i ^ + t r = -o x - q U 1 + iV)(x - q + ip ■ x, (2.14) 

1 — a 2 I — a 2 
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we can compare our results with those of Arick et al j3j. Remembering that 
a(t) = ae~ 2lt and ft(t) = fte~ %t (taking a and ft real in the rest of this section 
without loss of generality), it follows that the position q and momentum p 
of the wave-packet center are: 

rftcost ftsint^ 

q = 17^— > —i }> 2 - 15 

1 1 + a 1 — a } 



■ftsint/ 2a cos 2 1 \ ft cos t ^ 2a sin 2 1 
•TTa' 
Similarly, we find 



r /j sin t / z,a cus t \ cus l ( zu Sill b \ i 

P = tTT^i 1 + |l- tt e- 2 ^l 2 ) ' T^i 1 - \\-ae-^ )l- ^ 



«e 21 * 12 



and 



U = J- ^1 (2.17) 

1 — or 

2a sin 2t 

these matrices determine the probability density, 

p(x) = exp{ — (x — q)U _1 (x — q)} 
and the current density 

j(x) = p(x)[p + V(x- q)] 

(see Ref. 0). 

The above relations describe a Gaussian wave-packet with a periodically 
varying variance, and a center moving along the classical trajectory of a 
harmonic oscillator. This trajectory is an ellipse of eccentricity 



a. 



1 + a 

and with semi-major and semi-minor axis determined by the prameter ft : 

ft 

1 ±a ' 

Following the formulas of Ref. we can readily calculate the dispersions: 



which show that the width of the Gaussian wave-packet oscillates sinusoidally 
between the values |(1 — a)/(l + a) and ~(1 + a)/(l — a). Similarly, 



, .,, .,, ... 2a 2 sin 2 2t 1 1 . . 

= = ; ^m-^p + 2 u-^i. ' (2 - 20) 



a 



and 

, ,. ,. v , ,. ,. v Q! sin 2t , . 

(AxA Px .) = (AyA Pj/ ) = y— (2.21) 

From these formulas, we notice that: 

(Ax 2 )(Apl) - (AxAp x ) 2 = (Ay 2 )(Apl) - (AyAp y ) 2 = ~ , (2.22) 
which is precisely the definition of a coherent correlated state |2j. 



3 Polar operators and coherent states 

In order to compare with the standard results of the two-dimensional har- 
monic oscillator, we consider the usual creation and annihilation operators 
in Cartesian coordinates: 

i . d , 

at-^fr-A), (3.23) 

with (xi,x 2 ) = (x, y) = (r cos0, rsin0). They satisfy the commutation 
relations: 

Define now the operators 

A=^(a! + m 2 ) (3.24) 
v2 

B = -^(oi-ioz). (3.25) 
v2 



Their explicit forms in polar coordinates are: 

d_ 

dr r 



o V Br r rich 
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(3.26) 



2 v dr r i 
and they satisfy the commutation relations: 

[A, B]=0= [A, fit] 

[A, A*] = 1 = [B, St]. 

If we now express the wave functions in coordinate representations as in 
Eq. (|2.13J) . then it is easy to see that: 

1 + a 



(3.27) 



\f)r 



1 — a 



(3.28) 



and 



d id 
dr r dd 



1 + a 
1 — a 



a 



(3.29) 



These last two relations can be written in a coordinate independent form 
in terms of the operators defined in Eqs. ()3.26|) . Thus, for the righthanded 
polarized states: 



(A + aB^)\a(3)+ = 
(aA* + B)\aP)+ = f3\a(3). 



(3.30) 
(3.31) 



From these expression it is easy to find the expectation values of the A 
and B operators in the righthanded states: 



■(A)- 



-(B)- 



a(5* 



1 — \a\ 

P 



1 



a 



^3.32) 



In a similar way, it can be seen that for the lefthanded polarized states: 
(aAt + S)|a/3}_ = (3.33) 



(A + aSt)|«/3). 



/3\a/3). 



(3.34) 



and the expectations values in these states are: 



■(A)- 



a 



ad* 
1 — \a\ z 

Notice that the eigenvalue equations for the |or/3)_ are the same for the \ct/3) + 
states if the interchange A <-> B is made. 

As for the expectation values of quadratic operators, some lengthy but 
straightforward calculations show that: 

+ (AU) + = J"™! , (3-36) 



(1 - \a 



+ (B ] B)+ = TT^TT^ , (3-37) 



1 — \a 



+{A ' B) + = - (i -ht ' (3 ' 38) 



a 2 /3* 2 

-(^ 2 ) + = n ^ , (3-39) 



1 — a 



-(ff 2 )+ = 7i (3-40) 



1 — a 



2\2 ' 
2 



-<^>+ = ~ n , • (3-41) 



(1 - |a 



The equivalent expressions for the expectation values in the left-handed 
state can be obtained directly by the interchange A <-> I? combined with 
+ (...>+--<•.•>- 

Defining the Hermitian operators Qa,b and -FU.b as = \{Qa,b + 

iPa,b)i h follows from the above formulas that 

±((AQ A<B ) 2 ) ± =± ((AP A)B ) 2 )± = 1. (3.42) 

Therefore, the states under consideration are indeed coherent states, that is, 
they have minimum dispersion with respect to the A and B operators. 

Finally, we recall that the operators and a\ also generate an U(2) 
algebra (see, e. g., [S]). In terms of the standard operators a« and a\ and the 
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new operators A and B, the components of an angular momentum operator 
J can be denned as: 

1 



J x + iJ y = ^[(4 + ia\)(a x + ia y )} 
2 ld d 2 Id 2 fid Id 2 



r dr dr 2 r 2 d<p 2 
i 
2 



-t 2 86 r d(f)dr 



J~ = — (a],a~ — ala, 



- (AAt-BBi) = -4r 
2 v ! 2 86 



and the number operator is 



N = a\a x + a\a y = A 1 A + B^B, 



(3.43) 



(3.44) 



(3.45) 



which is simply related to the Hamiltonian appearing in Eq. (|2.1jl : H 
N + 1. 



4 Results and conclusions 

We have performed a numerical calculations in order to illustrate the evo- 
lution of the states described in this paper. In Fig. 1, we present a state 
corresponding to left handed polarization with particular values of the pa- 
rameters a and /3. As expected, the Gaussian wave packet rotates along a 
ellipse and changes its dispersion periodically; the form of the wave-packet is 
recovered after each cycle, in accordance with the results obtained by Arickx 
et al. 'S\. A cross section of the (^^(r, 6)\ 2 is shown in the figure. 

An interesting feature of the states we have defined in this article is that 
they are coherent correlated states with respect to the cartesian position 
and momentum operators. However, a new set of position and momentum 
operators can be defined in polar coordinates, and these same states turn 
out to be linear combinations of the eigenstates of these operators, with the 
peculiarity of being truly coherent with respect to these operators, having 
the minimum dispersion allowed by the Heisenberg uncertainty principle. 

In conclusion, the use of polar coordinates permits to see some features 
that are not evident in cartesian coordinates. The relation between classical 
and quantum motions is made explicit including the rotational degrees of 
freedom that correspond to different polarizations. 



9 



Acknowledgments 

We are grateful to R. Jauregui for many fruitful comments and discussions. 

References 

[1] Glauber R J 1963 Phys. Rev. 131 2766. 

[2] Dodonov V V, Kurmyshev E V, and Man'ko V I 1980 Phys. Lett A 79 
150. 

[3] Arickx F, Broekhove J, and Van Leuven P 1989 Nuovo Cimento B 104; 
1991 Phys. Rev. A 43 1211. 

[4] Padgett M, Courtial J, and Allen L 2004 Phys. Today 57, n. 5, 35. 

[5] Frank A and Van Isacker P 1994 Algebraic Methods in Molecular and 
Nuclear Structure Physics, chap. 1, Wiley (New York). 

Figure caption 

Fig. 1- An horizontal cross section of the square modulus ^^(r, 0)| 2 with 
parameters a = 0.4 and /3 — 5, in two-dimensional coordinates space. The 
cut is at ^^(r, 0)| 2 = 0.05. The classical trajectory, an ellipse with eccen- 
tricity 2y/a/ (1 + a), is added to illustrate the corresponding classical motion. 
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